Dynamic compression-shear response of brittle materials with specimen recovery

A new configuration for compression-shear soft-recovery experiments is presented. This technique is used to investigate various failure mechanisms during dynamic multiaxial loading of an Al₂O₃/SiC nanocomposite and TiB₂. Velocity profiles of the target surface are measured with a variable sensitivity displacement interferometer, yielding normal and transverse velocity-time histories. A dynamic shear stress of approximately 280 MPa is obtained, in the Al₂O₃/SiC nanocomposite, for an imposed axial stress of about 3.45 GPa on a 540 m thick sample. This dynamic shear stress is well below the value predicted by elastic wave propagation theory. This could be the result of stress-induced damage and inelasticity in the bulk of the sample or inelasticity on the sample surface due to frictional sliding. To gain further insight into the possible failure mechanisms, an investigation of compression-shear recovery techniques, with simultaneous trapping of longitudinal and lateral release waves, is conducted.     

增韧环氧树脂的动态裂纹扩展研究

本文主要进行了环氧及增韧环氧树脂的断裂韧性及裂纹快速扩展的试验研究。试验过程中采用了ＧＬＣ－１型高速裂纹扩展测试仪来测试裂纹的扩展速度,得到在裂纹扩展过程中裂纹扩展速度曲线。本文结合不同的计算公式及有限元分析方法,讨论了各个确定断裂韧性公式的准确程度,发现传统的静态断裂韧性的分析方法所得到的结果偏大,有一定的危险性,建议使用试验与数值计算相结合的方法;同时还发现增韧不仅可以提高材料的静态和动态断裂性能,而且在裂纹扩展过程中可以起到减缓裂纹扩展的作用     

Seismic gyrotropy as a result of rock-microstructure dissymmetry

Seismic investigations show that rocks exhibit not only anisotropic properties but also gyrotropic ones. Recently, the phenomenological theory has been developed to describe the phenomenon ‘seismic gyrotropy'. The paper gives a possible explanation of the gyrotropic properties of rocks on a microlevel. A gyrotropic model of grainy rock with dissymetry in its microstructure is presented. The model is built on a principle ‘azimuthal turn plus translation'. The mathematical and physical formulations of the problem are discussed. An algorithm is constructed to find the displacements, their first and second derivatives with respect to spatial coordinates as well as stresses inside a grain, an element of the model. Mathematical modelling of elastic-wave propagation in the model considered enabled to demonstrate that the model possesses gyrotropic properties, i.e., it rotates on a polarization vector. Gyration constants are determined.     

销钉对FCP电位测试精度影响之BEM分析

本文用边界元法对销钉接触角及销钉孔半径对三种电流输入方式下紧凑拉伸试样的电位函数影响进行了分析。结果证明,人们在电位法测试中常用的电流输入方式容易引起误差,并非最佳方案.最佳方案是另一种同样方便的电流输入方式.     

阶梯压电层合梁的波动动力学特性

采用行波理论系统地研究了压电阶梯梁的自由振动分析以及强迫响应的分析方法. 基于分布 参数理论研究了压电阶梯梁的波传播特性，忽略柔性梁横向剪切和转动惯量的影响，给出了 梁的轴向和横向的简谐波解. 将压电阶梯梁离散化为单元，考虑压电片的刚度和质量的影响， 建立了节点散射模型. 应用位移连续和力平衡条件，推导了节点的波反射和波传递矩阵，在 此基础上，引入波循环矩阵的概念，给出波循环矩阵、波传递系数矩阵的确定方法. 应用波 循环矩阵可以有效地计算结构的固有频率. 另外，应用波传递系数研究了压电陶瓷作动器位 置对其驱动能力的影响. 得出两个主要结论：1)作动器靠近悬臂梁固定端将有较强的驱动 能力，悬臂梁边界反射行波产生弯曲消失波有利于增大压电波的模态传递系数；2)模态传 递系数与固有频率的灵敏度密切相关，波传递系数越大, 对应该处固有频率变化灵敏度越大. 另外，数值算例表明了行波方法比有限元方法具有更高的计算精度.     

Theoretical and experimental analysis of longitudinal wave propagation in cylindrical viscoelastic rods

Wave propagation in viscoelastic rods is encountered in many applications including studies of impact and fracture under high strain rates and characterization of the dynamic behavior of viscoelastic materials. For viscoelastic materials, both material and geometric dispersion are possible when the diameter of the rod is of the same order as the wavelength. In this work, we simplify the Pochhammer frequency equation for low and intermediate loss viscoelastic materials and formulate corrections for geometric dispersion for both the phase velocity and attenuation. The formulation is then experimentally verified with measurements of the phase velocity and attenuation in commercial polymethylmethacrylate rods that are 12 and in diameter. Without correcting for geometric dispersion, the usable frequency range for determining the phase velocity and attenuation for the rod is about , and about for the rod. Using the correction procedure developed here, it was possible to accurately determine the phase velocity and attenuation up to frequencies exceeding for the rod and for the rod. These corrections are applicable to many polymers and other viscoelastic materials. From thereon, the viscoelastic properties of the material can be determined over a wide range of frequencies.     

Application of scaled boundary finite element method in static and dynamic fracture problems

The scaled boundary finite element method (SBFEM) is a recently developed numerical method combining advantages of both finite element methods (FEM) and boundary element methods (BEM) and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors (SIFs) directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.The project supported by the National Natural Science Foundation of China (50579081) and the Australian Research Council (DP0452681)The English text was polished by Keren Wang.     

VISCOPLASTIC SOLUTION TO FIELD AT STEADILY PROPAGATING CRACK TIP IN LINEAR- HARDENING MATERIALS

An elastic-viscoplastic constitutive model was adopted to analyze asymptotically the tip-field of moving crack in linear-hardening materials under plane strain condition. Under the assumption that the artificial viscosity coefficient was in inverse proportion to power law of the rate of effective plastic strain, it is obtained that stress and strain both possess power law singularity and the singularity exponent is uniquely determined by the power law exponent of the rate of effective plastic strain. Variations of zoning structure according to each material parameter were discussed by means of numerical computation for the tip-field of mode Ⅱ dynamic propagating crack, which show that the structure of crack tip field is dominated by hardening coefficient rather than viscosity coefficient. The secondary plastic zone can be ignored for weak hardening materials while the secondary plastic zone and the secondary elastic zone both have important influence on crack tip field for strong hardening materials. The dynamic solution approaches to the corresponding quasi-static solution when the crack moving speed goes to zero, and further approaches to the HR (Hui-Riedel) solution when the hardening coefficient is equal to zero.     

Two-dimensional transient dynamic response of orthotropic layered media

In this study, two-dimensional transient dynamic response of orthotropic plane layered media is investigated. The plane multilayered media consist of N different generally orthotropic, homogeneous and linearly elastic layers with different ply angles. In the generally orthotropic layer, representing a ply reinforced by unidirectional fibers with an arbitrary orientation angle, the principal material directions do not coincide with body coordinate axes. The solution is obtained by employing a numerical technique which combines the use of Fourier transform with the method of characteristics. The numerical results are displayed in curves denoting the variations of stress and displacement components with time at different locations. These curves clearly reveal, in wave profiles, the scattering effects caused by the reflections and refractions of waves at the boundaries and at the interfaces of the layers, and also the effects of anisotropy caused by fiber orientation angle. The curves properly predict the sharp variations in the response at the neighborhood of the wave fronts, which shows the power of the numerical technique employed in the study. By suitably adjusting the elastic constants, the results for multilayered media with transversely isotropic layers, or layers with cubic symmetry, or isotropic layers can easily be obtained from the general formulation. Furthermore, solutions for some special cases, including Lamb’s problem for an elastic half-space, are obtained and compared with the available solutions in the literature and very good agreement is found. Preliminary version presented at the Second International Congress on Mechatronics (MECH2K3), Graz, Austria, July 14-17, 2003.